Momentum-dependent mean field based upon the 
Dirac-Brueckner approach for nuclear matter 

G. Q. Li and R. Machleidt 

Department of Physics, University of Idaho, Moscow, ID 83843, USA 

Abstract 

A momentum-dependent mean field potential, suitable for application in 
the transport-model description of nucleus-nucleus collisions, is derived in a 
microscopic way. The derivation is based upon the Bonn meson-exchange 
model for the nucleon-nucleon interaction and the Dirac-Brueckner approach 
for nuclear matter. The properties of the microscopic mean field are examined 
and compared with phenomenological parametrizations which are commonly 
used in transport-model calculations. 
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The microscopic description of nuclear matter, finite nuclei, and nuclear reactions in 
terms of the realistic nucleon-nucleon (NN) interaction continues to be an interesting topic 
in nuclear physics. Although quantum chromodynamics (QCD) is believed to be the ultimate 
theory of strong interaction, the only quantitative NN potential models available up till now 
are based on the idea of meson exchange; a well-known example is the Bonn potential 

Special many-body theories, such as the Brueckner approach and the variational method, 
have been developed such that a realistic NN interaction can be applied in nuclear many- 
body calculations. Nuclear matter saturation has been quantitatively explained by the Dirac- 
Brueckner-Hartree-Fock (DBHF) approach starting from a realistic NN interaction. 

Thus, this approach provides a natural starting point for the self-consistent description 
of nuclear matter, nuclear structure, and nuclear reactions in terms of the realistic NN 
interaction. 

The extension of the DBHF approach from nuclear matter to the structure of finite 
nuclei and nucleon-nucleus scattering have been attempted 0-0]- The direct solution of 
the DBHF equation in finite nuclei is, however, rather involved. More practically, one either 
parametrizes the DBHF results in nuclear matter in terms of simple Lagrangians [§-0 
that can be easily applied in finite nuclei, or one performes a full Brueckner-Hartree-Fock 
(BHF) calculation in the finite nucleus taking the relativistic effects via the local-density 
approximation into account 

In tune with the latter approach, we present in this paper a momentum-dependent 
mean field single-particle potential derived from DBHF nuclear matter calculations. This 
mean field potential is suitable for application in the transport-model description of nucleus- 
nucleus collisions. Together with the in-medium NN cross sections derived in our earlier work 
13| , a fully self-consistent calculation of nucleus-nucleus collisions can then be performed. 
A self-consistent approach has been pursued in Refs. fT^-Q applying the Reid soft- 



core potential in a non-relativistic Brueckner-Hartree-Fock (BHF) calculation. The 
major criticism []T^ of this study comes from the fact that the BHF calculation cannot 



even reproduce quantitatively the saturation properties of static nuclear matter pO[; the 



application of its results to colliding nuclear matter is thus debatable. 

From the theoretical point of view, nucleus-nucleus collisions at intermediate energies 
are of particular interest. At these energies, both the mean field and two-body collisions 
play an equally important role in the dynamical evolution of the two colliding systems; 
therefore, they have to be taken into accout on an equal footing in transport models such 
as the Boltzmann-Uehling-Uhlenbeck (BUU) equation |T9| or quantum molecular dynamics 
(QMD) p. 

In the early applications of transport models to nucleus-nucleus collisions, the mean field 
and the in-medium NN cross sections were parametrized separately: a simple Skyrme-type 
parametrization was used for the mean field [0 , while for the in- medium NN cross sections 
free-space proton-proton (pp) data were applied p2| , p!9| . By varying the mean field (i.e., the 



parameters in the Skyrme parametrization) — which is related to the nuclear equation of 
state — and by comparing theoretical results with experimental data, one expects to obtain 
some information on the nuclear equation of state. 

There are at least two uncertainties in this approach: the mean field and the in-medium 
NN collisions. Some observables that are believed to carry useful information about the 
nuclear equation of state are also strongly affected by what is assumed for the NN collisions. 
Therefore, the in-medium NN cross sections used in transport models should be determined 
consistently before reliable information concerning the equation of state can be obtained. 
Within the DBHF approach, the in-medium NN cross sections below pion threshold have 
been studied in Ref . [|1^ . Using those microscopic in- medium NN cross sections will reduce 
the uncertainty concerning the in-medium NN collisions. 

Furthermore, it is well known that — due to correlations and exchange terms — the 
effective nuclear interaction is momentum dependent. Empirically, the momentum depen- 
dence of the mean field is best illustrated by the nucleon effective mass m*/m which is about 



0.7 at low energies and approaches unity at high energies Momentum dependence of 
the mean field has been found to play an important role in observables that are claimed to 
carry useful information about the nuclear equation of state, such as transverse momentum 



transfer p^j25[] and particle production ||2^ in nucleus-nucleus coUsions. The microscopi- 



cally derived momentum dependence of the mean field may help reducing also this source 
of uncertainty. 

Our derivation of the mean-field potential is based upon the DBHF approach for nuclear 
matter. For a comprehensive and pedagogical introduction into the basic ideas of the Dirac- 
Brueckner approach and its formalism, we refer the interested reader to section 10.5 of 
Ref. 1^ and Ref. Here, we will only briefly summarize the major points and the basic 
equations. For the NN potential, we will apply the one-boson-exchange representation of the 
Bonn meson-exchange model for the NN interaction |^ . It includes six nonstrange mesons 
with given mass and coupling. Pseudovector (derivative) coupling is used for pseudoscalar 
mesons (vr and rj). A form factor of monopole type is applied at each meson- nucleon vertex 
which simulates the short-range dynamics of quark-gluon nature. The potential provides an 
accurate description of the deuteron and NN scattering [@,|]. 

The essential point of the DBHF approach is the use of the Dirac equation for the 
description of single-particle motion in the nuclear medium 

[a-p + (3{m + Us) + Uv]u{p, s) = eu{p, s) (1) 

where Us is an attractive scalar field and Uy the time-like component of a vector field which 
is repulsive; m denotes the experimental mass of the free nucleon. 

The positive-energy Dirac spinors that solve Eq. (1) can be written as 



h( 1 ^ 



Ep + m^ ^ 



2m 



Xs (2) 



p 

\ Ep+rh I 



with 



fh = m + Us , (3) 
^p = (m' + p')^ (4) 



and Xs a Pauli spinor. The normalization is uu 
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Similar to conventional Brueckner theory, the basic quantity in the DBHF approach is 
a G-matrix which satisfies the in-medium Thompson equation (relativistic Bethe-Goldstone 
equation) [gg^ 



G(q', q; P) = V(q', q) + P / ^,V{c^, 



^2 



m 



(27r)3 ^Ef, 



{l/2)P+k 

^^^^1^^ GiKc-P) (5) 

For a pair of interacting nucleons with momentum pi and p2 in the nuclear matter rest 
frame, the center-of-mass (cm.) momentum is P = pi + p2 and their relative momentum 
q = (1/2) (pi — P2). V denotes the principal value. Note that Eq. (5) is density dependent 
due to the Pauli projector, Q, and the scalar field. Us; in the notation used in Eq. (5), this 
density dependence is suppressed. Notice also that V is density dependent due to the in- 
medium Dirac spinors, Eq. (2), representing the four outer legs of the one- meson-exchange 
Feynman diagrams defining the potential. The in-medium Dirac spinors depend on m which 
through Us depends on the density of the medium. 

The Dirac equation and the in-medium Thompson equation are solved self-consistently. 
The nuclear matter properties are derived from the G-matrix. For more details, see section 
10.5 of Ref. 11 and Ref. 

According to Eq (1), the single-particle energy of a nucleon with momentum pj is 

ei = ^ < i\j-Pi + m|i > +^Us + Uy (6) 
Ei Ei 

= ^ + ^Us + Uv (7) 

= E, + Uv (8) 

In our notation, \i > is a single particle state represented by a Dirac spinor of the kind 
Eq. (2) and < i\ is the adjoint Dirac spinor, u = u^jo', < i\i >= 1; = ct] P = 7o- 
The scalar and vector field. Us and Uy, are determined from 

Us + Uv = ReY. ^< uIGKj - 3^ > (9) 
Ei j<kp Eitjj 



This is the relativistic analogue to the non-relativistic Brueckner-Hartree-Fock definition of 
a single-particle potential using the effective mass approximation. Self-consistency of the 
DBHF calculations means satisfying Eq. (6). 



As shown in Ref. |3^, it is a good approximation to assume Us and Uy to be momentum 
independent (note, however, that they are density dependent). Notice also that momentum- 
independence of Us and Uy does not imply momentum- indepence of the mean field (cf. 
Eq. (13), below). Momentum-dependence of Us and Uy is one possible source of momentum- 
dependence of the mean field, but it is the less important one (see discussion below). 

In order to define a mean- field potential, we rearrange the single-particle energy, Eq. (8), 

ei = Ei + Uv + E,-Ei = Ei + U (10) 

with the free energy 

E, = {p^ + mr^\ (11) 

and the mean-field potential 

U, = Uv + E,-Ei. (12) 

More explicitly, this mean field is 

U{p,p) = Uvip) + [{m + Us{p)f + pT^' - + V'f" (13) 

where the dependence on the density p is now clearly indicated (p = |p|). Us{p) and Uv(p) 
are listed for many densities in Table VII of Ref. 0; e. g., for po = 0.17 fm~^ one has 
Us = —355.7 MeV and Uy = 274.7 MeV. By using this table in a spline interpolation 
(adding Us = Uy = for p = 0), one can obtain Us and Uy with good accuracy for any 
density between and 4.2po- For a simple parametrization of our mean field potential, see 
Eq. (20), below. 

The mean field, Eq. (13), shows strong momentum dependence for small p. For large 
p {p » m), the square roots are ~ p and, thus, they cancel. Therefore, U vanishes 
approximately for large p (notice that Uy can be neglected when p becomes very large). In 
summary, U behaves in a physically reasonable way. 



As mentioned, we have neglected the momentum dependence of Us and Uy- In Ref. [^] it 



is shown that, for low momenta, the momentum dependence oi Us and Uy can be described 
by 

f/.(p)=Cr-C,?'^ (14) 

Pf 

Uy{p) = uP -U^^^ (15) 
Pf 



with pp the Fermi momentum. It is found that P,|30 



^ - ^ -0.05 (16) 



Expanding Eq. (13) up to terms in and using Eqs. (14)-(16), one obtains 



r..--0.05^^i^/.!|=|^p^ ,17, 



The last term on the r.h.s. describes the momentum dependence of Eq. (13) for low mo- 
menta; the second but last term (which is typically a factor of five smaller) is the additional 
momentum dependence from Us and Uy that is obviously negligible. Notice also that Us 
and Uy are of opposite sign, but of the same order of magnitude. 

The nucleon effective mass, m*, as predicted by the DBHF calculations for nuclear matter 
is discussed in detail in Ref. . Using the Bonn A potential, which we apply in this work, 
one gets m* /m = 0.71 at normal nuclear matter density (po = 0.17 fm^'^) — a reasonable 
value. 

Empirically, the most sophisticated investigation of the nucleon mean field (optical po- 



tential) is conducted in Dirac phenomenology p4| , p6| . In this approach, the scalar and vector 



potentials are parametrized with a number of free parameters fitted to nucleon-nucleus scat- 
tering data. In order to compare our microscopic mean field obtained in DBHF calculations 
with the empirical one from Dirac phenomenology, we extract from the Dirac global code 
p6| the scalar and vector potentials in the centers of ^°Ca and Pb. In Fig. 1 we com- 
pare our microscopic mean field (solid curve) at normal nuclear matter (po=0.17 fm^^) with 
those extracted from the Dirac phenomenology. The open circles and squares correspond to 



empirical values in the center of ^o^Pb based on fit 1 and fit 3 of Ref. ||36|, respectively, while 



the solid triangles are the empirical values in the center of ^°Ca based on fit 1 of Ref. |^ 
The difference between the circles and squares reflects the uncertainties in the empirical de- 
termination of the mean fleld, while the difference between the circles and triangles shows, 
to some extent, flnite-size effects. In general, our microscopic result is in good agreement 
with the empirical one. The microscopic mean fleld, shows, however, a stronger momentum 
dependence. This may be due to the assumption of momentum independence of the scalar 
and vector flelds in our DBHF calculation. Finite-size and surface effects, included in the 
Dirac phenomenology, could also play a role. 

In the past, phenomenological parametrizations of momentum-dependent mean- fleld po- 
tentials have been constructed and used in transport models. We mention here two typical 



examples. One is proposed in Ref. p5[ and used in BUU calculations [|T^, we denote it 



by U^^^; the other is proposed in Ref. and used in QMD calculations I^TJ, we denote 



it by For U^^\ one parameter set has been given that corresponds to a soft 

equation of state with incompressibility K = 215 MeV; we denote this parameter set by 
GBG. For U'^'^\ two parameter sets, corresponding to a soft and hard equation of state, have 
been proposed 0,^; they are usually denoted by SMD {K=200 MeV) and HMD {K=380 



MeV), respectively. In static nuclear matter, U^^^ and f/*-^-* are given by 

U^'\p,p) =ap + 13 p' + 5(0.7965 + Y^^TfJ^^P ^^^^ 

U^^\p,p) = ap + f3p^ + (51n2(£p2 + l)p- (19) 

Parameter sets for f/(^) (GBG) and f/^^^ (SMD, HMD) are hsted in Table 1. 

In Fig. 2, we show the density dependence of our microscopic (Eq. (13)) and of phe- 
nomenological (Eqs. (18) and (19)) mean flelds. At low densities, the four mean flelds shown 
are very close to each other. At high densities, there are signiflcant differences: the HMD 
mean fleld is much stiffer than the others; the microscopic mean fleld is more close to the 
soft (SMD and GBG) phenomenological mean flelds. There is also an important difference 
between the microscopic mean fleld obtained in our (relativistic) DBHF calculation and the 
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one from non-relativistic BHF: the former is stiffer than the latter (see Fig. 4 of Ref. |T5 
and Fig. 3 of Ref. |]1^) due to repulsive relativistic effects. 

In Fig. 3, we show the momentum dependence of the microscopic and phenomenological 
mean fields. Both microscopic and phenomenological mean fields increase with increasing 
momenta. For the momenta shown, the microscopic mean field has a stronger momentum 
dependence than the phenomenological ones. Such a strong momentum dependence is also 
observed in the microscopic mean field from non-relativistic BHF calculations Note, 



however, that for very large momenta (which are not shown in Fig. 3) the situation is re- 
versed: the momentum- dependence disappears from our microscopic mean field (as discussed 
above), while the phenomenological parametrization continues to grow. 

For convenience in applications, we have parametrized our microscopically derived mean 
field using exactly the same ansatz as for f/*^^\ Eq. (19), namely, 

U{p,p)=ap + pp^ + 6\n\ep^ + l)p'' (20) 

The parameters are listed in Table 1, row 'DBHF'. The quality of the fit can be seen in Fig. 
4 where solid curves are obtained from Eq. (20) while the circles correspond to the exact 
calculation using Eq. (13). Both density dependence (upper part of Fig. 4) and momentum 
dependence (lower part of Fig. 4) are reproduced well by the parametrization Eq. (20), for 
the densities and momenta shown. Note, however, that the parametrization Eq. (20) is bad 
for large p {p >> m) since it continues to grow with increasing p, while the exact mean 
field, Eq. (13), becomes independent of p for large momenta. These very large momenta 
may, however, not be important in nucleus-nucleus collisions at intermediate energies. 

The incorporation of the microscopic mean field and the in-medium NN cross sections in 
transport-model calculations is straightforward. The propagation of nucleons is determined 



by Hamilton's equations of motion |[T9|j21 



^ = -^'U (21a) 



where the position dependence of the mean field enters through the position dependence 
of the density via the local-density approximation. Moreover, using microscopic in-medium 



NN cross sections ||T3[, the in- medium collisions can be treated in the same way as in usual 
transport models. The self-consistent investigation of nucleus-nucleus collisions based upon 
realisitic NN interactions will be reported elsewhere. 

In summary, we have derived a microscopic momentum dependent mean field using the 
DBHF approach for nuclear matter and the Bonn potential. The momentum dependence of 
our mean field comes out very close to that obtained in Dirac phenomenology. Qualitatively, 
our microscopic mean field shows similar density and momentum dependence as some of the 
phenomenological parametrizations, which are often used in BUU and QMD calculations; 
quantitatively, however, there are important differences, especially at high densities and 
momenta. For practical purposes, we have parametrized our microscopic mean field in terms 
of a simple analytic form which accurately reproduces our exact results; this can easily be 
applied in transport-model calculations. 
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Table 1: Parameters of phenomenological and microscopic mean fields. The parameter 
set 'DBHF' is the parametrization of our microscopic mean field, Eq. (13), in terms of the 
ansatzEq. (20). GBG uses Eq. (18), while SMD and HMD use Eq. (19). 



a 


(MeV-fm^) p (MeV-fm^T) 


7 


S (MeV-fm3<^) 


e (MeV-2) 


a 


A (MeV) 


GBG 


-888.96 


1687.53 


7/6 


-460.12 






400 


SMD 


-2294.12 


2412.34 


1.14 


9.24 


5.0x10-^ 


1.0 




HMD 


-764.71 


2394.49 


2.09 


9.24 


5.0x10"^ 


1.0 




DBHF 


-1129.51 


1501.31 


1.48 


35.07 


4.3x10-^ 


0.7 
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FIGURES 

FIG. 1. Momentum dependence of empirical and microscopically derived mean field potentials. 

The microscopic mean field is calculated for nuclear matter with p = 0.17 fm^'"^. The empirical 
mean fields are extracted from Dirac phenomenology at the centers of ^°Ca and ^°®Pb, the open 
circles and squares correspond to ^osp]-, based on fit 1 and fit 3 of Ref. [36], while the solid triangles 
correspond to ^''Ca based on fit 1. 

FIG. 2. Density dependence of microscopic and phenomenological mean fields for several nu- 
cleon momenta. The solid curve represents our microscopic mean field as obtained from DBHF. 
The dashed and dash-dotted lines are phenomenological parametrizations proposed in Ref. [24]. 
The dotted line is from Ref. [19,25]. 

FIG. 3. Momentum dependence of microscopic and phenomenological mean fields for several 
densities. Notation as in Fig. 2. 

FIG. 4. Comparison of the parametrization Eq. (20) (solid curves) with the exact results (dots) 
for our microscopic mean field, Eq. (13). The upper part shows the density dependence, and the 
lower part the momentum dependence. 
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